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Solar  thermoelectric  generators  (STG)  using  cheap  parabolic  concentrators 
with  high-ZT  modules  can  be  a  cost-effective  alternative  to  solar  photovoltaics 
for  micropower  generation.  A  thermodynamic  analysis  is  presented  for 
predicting  the  thermal-to-electrical  conversion  efficiency  for  the  generator. 
With  solar  concentration  of  66  x  suns,  a  system  efficiency  of  3%  was  measured 
for  a  commercial  Bi2Te3  module  with  output  power  of  1.8  W.  Using  novel 
thermoelectric  materials  such  as  ft-type  ErAs:(InGaAs)i_JC(InAlAs)JC  and p-type 
(AgSbTe^PbSnTe)!.*,  a  conversion  efficiency  of  5.6%  can  be  achieved  for  a 
STG  at  120  x  suns. 
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INTRODUCTION 

Historically,  thermoelectrics  have  been  used  pri¬ 
marily  for  deep-space  exploration  and  waste  heat 
recovery.  We  explore  the  potential  of  thermoelec¬ 
trics  with  solar  energy  for  electricity  generation 
using  a  solar  thermoelectric  generator  (STG) 
(Fig.  la).  A  solar  collector  (parabolic  reflector) 
directs  the  sunlight  onto  a  fixed  focal  spot.  The  hot 
side  of  a  thermoelectric  generator  placed  at  this  spot 
heats  up  as  it  absorbs  the  concentrated  sunlight. 
A  “selective  surface”  on  the  hot  side  of  the  module 
has  high  (>90%)  absorbance  from  300  nm  to 
1000  nm.1  This  allows  the  hot  side  to  absorb  most  of 
the  incident  solar  energy.  With  appropriate  selec¬ 
tive  surfaces  most  of  the  solar  energy  [ultraviolet 
(UV),  visible,  and  infrared  (IR)]  incident  on  the 
generator  contributes  towards  raising  the  hot-side 
temperature.  A  fraction  of  this  thermal  energy  is 
converted  to  electrical  energy  by  the  thermoelectric 
module  via  the  Seebeck  effect.  This  paper  provides  a 
brief  overview  of  the  previous  work  in  this  field, 
followed  by  a  presentation  of  new  thermodynamic 
and  cost  models  for  STGs.  Finally,  experiments  are 
used  to  validate  the  thermodynamic  model;  these 
experiments  also  realize  record  performance  for 
STG  modules. 


(Received  July  10,  2009;  accepted  March  8,  2010; 
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The  first  concept  for  a  STG  was  tested  in  1922  by 
Coblentz  for  measuring  infrared  radiation  from 
stars.2  Early  work2-5  showed  low  system  efficiency 
(<  1%),  primarily  due  to  low  module  ZT  (<0.4),  low 
solar  concentration,  and  low  hot-side  temperature. 
To  date,  the  highest  measured  efficiency  for  solar 
thermoelectrics  is  3.35%,  using  a  unicouple  with 
ZT  =  0.4  and  solar  concentration  of  48  x  suns 
(1  sun  =  1000  W/m2).2  The  measured  output  power 
was  0.156  W.  Figure  2  summarizes  the  previously 
published  experimental  results  for  different  solar 
thermoelectrics,  along  with  the  experimental 
results  presented  at  the  end  of  this  paper. 

Complementing  this  experimental  effort  has  been 
a  history  of  STG  modeling.  A  recently  published 
model  estimates  a  maximum  system  efficiency  of 
35%  for  a  flat-plate  (lx  solar  concentration)  solar 
thermoelectric  system.6  Herein  we  show  that 
the  recently  reported  models  are  unphysical — 
specifically  they  treat  the  hot-side  heat  transfer 
inconsistently,  resulting  in  hot-side  temperatures  of 
900  K  without  solar  concentration.6  In  fact,  for  a 
Bi2Te3  module  with  natural  air  convection,  the  heat 
flux  at  low  solar  concentration  using  flat-plate  or 
box-type  concentrators  results  in  a  relatively  low 
hot-side  temperature  (<320  K). 

THERMODYNAMIC  MODEL 

We  present  a  thermodynamic  analysis  based  on 
energy  balance  and  heat  transfer  that  allows  us  to 
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Fig.  1.  (a)  A  solar  thermoelectric  generator,  showing  the  concentrator  and  the  heat  sink,  (b)  Schematic  of  the  enclosure  holding  the  thermo¬ 
electric  generator.  A  Fresnel  lens  is  used  to  increase  the  flux  concentration. 
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Fig.  2.  Previous  experiments  for  STG,  shown  with  system  efficiency 
and  output  electrical  power. 


predict  system  efficiency  ( rjsys )  for  the  STG.  The  loss 
mechanism  in  the  system  can  be  categorized  into 
two  parts:  loss  in  the  concentrator  and  loss  in  the 
thermoelectric  module.  The  concentrator  efficiency 
( rjQ )  is  defined  by  how  well  the  solar  flux  can  be 
guided  into  the  thermoelectric  module.  The  effi¬ 
ciency  is  given  by  r\ q  =  where  Qh  is  the  total 

thermal  power  into  the  thermoelectric  generator 
(TEG),  </>i  is  the  solar  flux  hitting  the  focal  spot  of  the 
concentrator,  and  A  is  the  area  of  the  generator.  The 
main  losses  degrading  the  concentrator  efficiency 
are  convective  and  radiative  losses.  The  optical 
absorbance  of  the  generator's  hot  side  is  also 
included  in  this  efficiency.  The  module  efficiency 
(t/teg)  depends  on  the  thermoelectric  material  and 
the  module  design.  It  is  given  by  rjTEG  =  T-,  where  P 
is  the  total  electrical  power  generated,  which 
depends  on  material  properties  such  as  the  Seebeck 
coefficient,  electrical  conductivity,  and  thermal 
conductivity.  For  maximum  power,  the  module 
efficiency  can  also  be  expressed  by  Eq.  1,  where  Th 
is  the  hot-side  temperature,  Tc  is  the  cold-side 
temperature,  and  AT  is  the  temperature  difference 
across  the  module.  ZT^^ \e  is  the  figure  of  merit, 
which  depends  on  the  material  parameters  (Seebeck 
coefficient,  thermal  conductivity,  and  electrical 
conductivity)  and  module  parasitics  as  described 
later  in  this  paper. 


AT 

'?TEG  “  2Th  -  AT/2  +  2(Th  +  Tc)/ZTmodule- 


(1) 


The  total  system  efficiency  is  given  by  rjsys  =  t/teg^q- 

In  general,  convective  loss  can  be  reduced  by 
introducing  suppression  mechanisms  such  as  mul¬ 
tiple  glass  panels  and  dead  air  between  the  glass 
and  the  hot  side.  In  our  experiment  (Fig.  lb),  the 
light  hits  the  generator  from  below.  An  enclosure 
traps  the  warm  air  near  the  hot  surface.  This 
trapped  air  suppresses  convective  loss  at  the  hot 
side.  For  radiative  loss  suppression,  the  hot  side  is 
coated  with  a  selective  surface  consisting  of  silicon 
polymer  as  a  binder  with  an  oxide  pigment.1  The 
selective  surface  has  a  large  absorbance  (0.88  to 
0.95)  for  visible  wavelengths  (300  nm  to  1200  nm) 
and  low  emissivity  (0.2  to  0.4)  at  wavelengths  above 
2.2  /im.1  The  selective  surface  allows  for  near- 
optimal  absorption  of  incident  solar  radiation  while 
minimizing  radiative  heat  loss.  Even  at  a  relatively 
low  hot-side  temperature  (<500  K),  the  selective 
surface  performs  better  than  ordinary  black  paint 
by  increasing  the  overall  temperature  drop  across 
the  TEG  by  10%. 

STG  performance  was  modeled  using  mature 
thermoelectric  materials  such  as  micro-alloy  Bi2Te3 
and  SiGe,  as  well  as  novel  materials  with  high 
material  ZT  (>1.2)  such  as  ErAs:(InGaAs)i_*- 
(InAlAs)*7  and  (AgSbTe)JC(PbSnTe)i_Jc.8  Thermo¬ 
electric  modules  have  relatively  low  effective  ZT 
(0.25  to  0.4)  compared  with  material  ZT  (Bi2Te3: 
ZT  «  1),  mostly  due  to  thermal  and  electrical 
parasitics.  Material  ZT  is  given  by  ZTm ateriai  = 

(#te^te)  ^  where  a  is  the  Seebeck  coefficient,  i?TE 

is  the  electrical  resistance,  and  KTE  is  the 
thermal  conductance  for  a  thermoelectric  material. 
Module  ZT  can  be  formulated  as  ZTm0(^u \e  = 

(  Tr  Tp - \(K  +JT - y)T’  where  ^parasitic  and 

^parasitic  are  the  electrical  parasitics  (primarily 
metal  interconnects  between  couples)  and  thermal 
parasitics  (primarily  ceramic  plates  for  conventional 
designs).  It  can  be  seen  from  the  above  relationships 
that  the  module  ZT  is  less  than  the  material  ZT  due 
to  parasitics.  The  module  conversion  efficiency  is 
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Fig.  3.  Depiction  of  heat  flux  through  the  generator,  showing  various 
interfaces  and  temperatures. 


determined  by  the  module  ZT  rather  than  the 
material  ZT. 

For  the  simulation,  the  set  of  Eqs.  2-9  repre¬ 
senting  heat  transfer  and  energy  balance  within 
the  STG  system  (Fig.  3)  is  self-consistently  solved 
for  temperatures  at  various  interfaces.  The  input 
solar  flux  (</>i)  is  the  only  independent  variable  in 
the  simulation.  The  cold-side  heat  transfer  coeffi¬ 
cient  ih)  is  set  by  the  natural  air  convection, 
which  limits  the  cold-side  temperature  for  the 
generator  (Tc).  The  thermal  impedance  of  the 
module  sets  the  hot-side  temperature  for  different 
heat  flux.  Depending  on  the  temperature  of  the 
module  (both  hot  and  cold  side),  the  generator 
efficiency  is  calculated.  Both  the  convective  and 
radiative  losses  at  the  hot  side  depend  on  the  hot- 
side  temperature. 

Qh  =  ~  ki(Th  -  Ta)A  -  k2T^A,  (2) 

where  rjQ  is  the  optical  absorbance  for  the  hot-side 
absorber,  k1  is  the  convective  heat  transfer  coeffi¬ 
cient,  k2  is  the  radiative  heat  transfer  coefficient, 
and  Ta  is  the  ambient  temperature. 

T1  and  T2  are  the  temperatures  across  the 
material.  Variations  in  these  temperatures  from 
Th  and  Tc  are  due  to  the  thermal  parasitics,  which 
are  captured  in  the  model  as  &h  and  kc — the 
hot-side  and  cold-side  thermal  conductance 
(-0.67  W/m2  K). 


Qh  =  MThA  Tl).  (3) 

The  total  heat  transported  from  the  hot  side  can 
be  separated  into  contributions  from  the  Seebeck 
effect,  thermal  conductance,  and  Joule  heating. 


a  is  the  Seebeck  coefficient  of  the  TE  material, 
which  is  temperature  dependent.  K  is  the  thermal 
conductance,  and  R  is  the  resistance  of  the  module. 
I  is  the  electrical  current  through  the  TEG. 

I2R 

Qh  =  a  TV  +  K(T±  -  T2)  -  — .  (4) 

Similarly,  the  total  heat  from  the  cold  side  is  also 
due  to  the  Seebeck  effect,  thermal  conductance,  and 
Joule  heating: 

T2R 

Qc  =  olT2I  +  JT(Ti  -  T2)  +  — .  (5) 

The  total  electrical  power  generated  by  the  TEG  is 
then 


p  =  Qh  -  Qc  =  <T  1  -  T2)I  -  i2r  =  I2Rl.  (6) 

For  maximum  power,  the  condition  of  load  matching 
is  applied,  according  to  which  the  load  resistance 
CRl)  is  equal  to  the  resistance  of  the  TEG  (R),  i.e., 


The  heat  transfer  at  the  cold  side  is  given  by  Eqs.  8 
and  9.  Ahs  is  the  area  of  the  heat  sink. 


n  kc(T2-Tc) 
Qc  = - t - , 


(8) 


Qc  =  hA\^{Tc  -  Ta).  (9) 


Using  this  analysis,  Bi2Te3  with  module  ZT  of 
0.64  gives  a  system  efficiency  of  4%  using  natural 
air  convection  at  the  cold  side  with  solar  flux  of 
70  x  suns.  In  the  thermodynamic  analysis,  the 
material  properties  (thermal  conductivity,  electrical 
resistivity,  and  Seebeck  coefficient)  are  taken  as 
temperature-dependent  parameters.  Figure  4 
shows  simulations  of  module  ZT  and  system  effi¬ 
ciency  at  various  input  solar  fluxes.  In  this  result, 
the  existence  of  peak  efficiency  is  due  to  the  trade¬ 
off  between  the  hot-side  temperature  and  the 
material  performance  at  high  temperature:  from 
Eq.  1,  the  module  efficiency  increases  as  the  hot-side 
temperature  rises.  However,  as  the  temperature 
increases,  Bi2Te3’s  ZT  degrades.  The  temperature  of 
900  K  in  Ref.  6  is  not  achievable  primarily  because 
the  radiative  loss  (for  s  =  1)  exceeds  the  incident 
solar  flux  and  because  the  material  ZT  for  Bi2Te3 
degrades  at  such  high  temperature. 

EXPERIMENTAL  PROCEDURES 

Experiments  were  conducted  using  a  Newport 
solar  simulator  (2  inch  x  2  inch  beam  size)  with 
commercial  TEG  modules.  The  average  input 
thermal  flux  into  the  module  was  66  x  suns,  which 
was  measured  by  using  a  thermopile  sensor.  This  is 
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Solar  flux  (  x  suns)  [1  sun  =  1 000W/m2] 

Fig.  4.  Theoretical  simulation  for  Bi2Te3,  showing  module  ZJ  and 
system  efficiency;  without  any  parasitics,  a  maximum  system  effi¬ 
ciency  of  4%  is  achieved  at  70 x  suns.  The  model  is  verified  with  an 
experiment  at  66x  suns.  Parasitics  are  included  in  the  theory  to 
represent  the  commercial  module. 


the  input  flux  on  the  secondary  concentrator,  i.e., 
the  Fresnel  lens.  The  generator  (TG12-4)  from 
Marlow  Inc.,  with  a  module  ZT  «  0.419  and  mate¬ 
rial  ZT  «  0.9, 9  gave  a  system  efficiency  of  3%  with 
output  electrical  power  of  1.8  W.  The  generator  was 
connected  to  a  power  metal-oxide-semiconductor 
field-effect  transistor  (MOSFET)-based  variable 
resistance  for  load  matching  and  maximum  power 
measurement.  An  electrical  parasitic  of  1.5  Q  and 
thermal  resistance  of  0.35  K/W  for  the  ceramic 
plates  along  with  Bi2Te3  material  parameters  were 
used  to  model  the  commercial  module.  The  experi¬ 
mentally  measured  system  efficiency  for  the  STG 
agrees  with  the  thermodynamic  model  presented 
above.  However,  we  can  see  that,  even  with  greater 
solar  concentration  and  with  no  module  parasitics, 
the  system  efficiency  with  Bi2Te3  is  limited  to  4%. 
In  order  to  realize  higher  efficiencies  we  must 
consider  materials  with  high  ZT  at  higher  tem¬ 
peratures. 

Novel  thermoelectric  materials  such  as  72-type 
ErAs:(InGaAs)i_JC(InAlAs)JC7  and  p-type  (AgSbTe)*- 
(PbSnTe)!.*,8  have  shown  promising  high  material 
ZT  (>1.2)  for  average  temperatures  in  the  range 
from  700  K  to  750  K.  Using  this  novel  material 
combination,  a  system  efficiency  of  5.64%  can  be 
achieved  with  ZT  of  0.64  (Fig.  4)  at  an  average 
temperature  of  546  K  (Fig.  5).  The  solar  concentra¬ 
tion  needed  to  reach  an  average  temperature  of 
546  K  is  120  x  suns.  Simulations  were  not  per¬ 
formed  beyond  a  hot- side  temperature  of  700  K,  as 
the  materials  have  not  been  tested  at  these  higher 
temperatures.  Thus,  for  the  simulation,  ZT  never 
reaches  the  highest  value  of  1.2.  The  thermody¬ 
namic  model,  even  with  novel  materials,  suggests 
that  system  efficiencies  of  6%  are  beyond  the  reach 
of  STGs.  However,  even  with  these  modest  effi¬ 
ciencies,  STGs  may  be  of  commercial  interest  if  the 
overall  system  cost  can  be  kept  sufficiently  low. 


Solar  flux  (x  suns)  [1  sun  =  1000  W/m2] 

Fig.  5.  Theoretical  simulation  for  ErAs/LAST,  showing  system  effi¬ 
ciency  of  5.64%  with  module  ZJ  of  0.64.  The  hot-side  temperature  is 
700  K.  Simulation  was  not  done  beyond  this  temperature,  as  the 
materials  have  not  been  tested  at  higher  temperatures. 


COST  MODEL 

In  this  final  section,  we  present  a  cost  model  for 
the  STG,  comparing  it  with  photovoltaics  (PV). 
The  bulk  retail  price  for  a  15-W  thermoelectric 
generator  is  US  $25, 10  giving  a  peak  power  price  of 
US  $1.67AVP,  whereas  the  retail  peak  power  price 
for  a  PV  module  is  US  $4/Wp  at  the  time  of  publi¬ 
cation.11  This  PV  price  includes  installation  cost. 
Both  TEG  and  PV  have  similar  lifetime  (15  years  to 
20  years)  and  low  maintenance  cost.  The  electronic 
circuit  requirements  are  also  similar  for  PV  and 
TEG,  as  both  require  batteries,  charge  controllers, 
inverters,  etc.  The  critical  difference  is  the  conver¬ 
sion  efficiency,  which  is  between  7%  and  12%  for  a 
low-cost  PV  module.  Commercial  TEG  modules  are 
rated  for  4%  optimum  conversion  efficiency  for  a 
large  temperature  difference  (TH  =  230°C, 
Tc  =  50°C).9  A  low-cost  parabolic  concentrator  can 
focus  a  solar  flux  of  ~10x  suns  on  a  spot  size  of 
28  cm  diameter.  Using  inexpensive  Fresnel  lenses 
as  secondary  concentrators,  the  flux  can  be  easily 
increased  to  50  x  to  150  x  suns.  The  electricity  price 
(US  0/kW  h)  for  the  STG  was  calculated  using  the 
total  amount  of  electrical  power  generated  within 
the  lifetime  of  a  generator  (20  years).  The  generator 
efficiency,  rj((p ±),  depends  on  the  incident  solar  flux 
(4> i)  as  shown  in  Fig.  6. 

The  output  power  at  any  given  instant  is  equal  to 
the  system  efficiency  multiplied  by  the  total  thermal 
power  incident  on  the  generator.  The  thermal  power 
into  the  generator  depends  on  the  incident  solar  flux 
and  the  generator  area  (A),  as  well  as  on  radiative 
and  convective  losses  from  the  absorber  surface. 
Figure  6  also  shows  the  total  number  of  hours  in  a 
year  (P^.)  that  solar  radiation  of  a  certain  intensity 
hits  the  Earth's  surface.  The  yearly  average  for 
this  particular  location  (Kathmandu,  Nepal)  is 
210  W/m2.  The  total  output  power  within  the  life¬ 
time  of  a  generator  is  given  by  the  sum  of  output 
power  for  each  hour  of  the  day  as  shown  in  Eq.  10. 
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Solar  flux  ( x  1000  W/m2) 

Fig.  6.  Dependence  of  system  efficiency  on  solar  flux  for  Bi2Te3  and 
ErAs/LAST  material.  The  total  number  of  hours  within  a  year  for 
certain  flux  radiation  is  also  shown  as  a  histogram. 


(^electrical)  =  (  ^  *?(&)  X  fa  X  A  X  P^.  j  X  20.  (10) 


Table  I  shows  the  electricity  price  with  the  STG 
using  various  thermoelectric  materials. 

We  explore  the  effect  of  parasitics  on  the 
electricity  pricing.  At  fixed  solar  concentration, 
the  module  ZT  is  varied  by  introducing  electrical 
parasitic  resistance  (0  Q  to  5  Q).  As  the  parasitic 
increases,  the  module  ZT  and  the  system  efficiency 
decrease.  When  the  system  efficiency  decreases,  the 
electricity  price  rises.  Figure  7a  shows  the  price 
dependence  for  Bi2Te3.  Simulations  with  two  dif¬ 
ferent  geometrical  aspect  ratios,  i.e.,  cross-section/ 
length,  of  0.56  cm  and  0.2  cm  for  the  TE  leg  show 
that  the  modified  Bi2Te3  module  with  low  para¬ 
sitics  and  longer  TE  legs  can  have  electricity  price 
comparable  to  the  PV  price.  For  higher  solar 
concentration  and  the  ErAs/lead-antimony-silver- 
tellurium  (LAST)  combination,  even  with  the  aspect 
ratio  of  0.56  cm,  which  is  the  aspect  ratio  for  a 
commercial  module,  the  system  efficiency  is  high 
enough  for  the  electricity  cost  to  be  within  the  PV 
price  range  (Fig.  7b). 


CONCLUSIONS 

In  this  work,  we  have  presented  experimental  and 
theoretical  results  relating  to  solar-to-electrical 
conversion  using  thermoelectrics.  Generators 


(a) 


Fig.  7.  (a)  Electricity  cost  variation  in  Bi2Te3  by  changing  module 
parasitic  resistance  at  70 x  suns.  As  the  resistance  increases,  the 
module  ZT  decreases.  With  a  modified  module  design  (longer  TE 
legs),  the  system  efficiency  can  be  increased  to  decrease  the  elec¬ 
tricity  cost,  (b)  By  choosing  higher  solar  concentration  (120x  suns) 
and  materials  such  as  ErAs/LAST,  the  electricity  cost  is  within 
22  US  0/kW  h  to  30  US  0/kW  h. 


employing  light  concentrators  and  modules  with 
high-ZT  thermoelectric  materials  are  an  attractive 
alternative  to  solar  photovoltaics  for  micropower 
applications.  We  can  predict  conversion  efficiency 
for  a  STG  using  various  materials  at  different 
temperature  ranges  using  a  thermodynamic  analy¬ 
sis  based  on  energy  balance  and  heat  transfer.  With 
novel  thermoelectric  materials  which  have  high 
material  ZT,  a  system  efficiency  of  5.6%  can  be 
achieved.  Experimentally  we  show  a  system  effi¬ 
ciency  of  3%  with  a  commercial  Bi2Te3  module, 
which  is  (to  the  best  of  our  knowledge)  the  highest 


Table  I.  Price  of  electricity  for  a  STG  using  various  thermoelectric  materials  compared  with  PV 


Material 

ZT 

(module) 

Max.  System 
Efficiency  (%) 

Solar 

Concentration 

Average 

Temperature  (K) 

Electricity  Price 
(US  0/kW  h) 

Bi2Te3 

0.64 

4 

70  x  suns 

446 

35 

ErAs/AgSbTe 

0.64 

5.64 

120  x  suns 

546 

22 

Si  (PV) 

- 

11 

- 

- 

20-2611 

1740 


Amatya  and  Ram 


efficiency  for  solar  thermoelectrics  using  a  module. 
From  the  cost  analysis  perspective,  even  though 
present-day  commercial  modules  have  low  module 
ZT  (^0.4)  and  low  system  efficiency,  the  thermo¬ 
electric  element  geometry  can  be  optimized  for 
performance  in  the  STG  application.  With  better 
module  design  (longer  TE  legs  and  low  parasitics), 
and  by  using  materials  with  improved  thermoelec¬ 
tric  properties  at  higher  temperature,  solar  ther¬ 
moelectrics  can  be  economically  competitive  with 
small-scale  PV  power  generation  for  many  applica¬ 
tions  such  as  rural  electrification  in  developing 
countries  and  power  supply  for  remote  sensors. 

REFERENCES 

1.  SOLKOTE:  HI/SORB-II  Selective  Solar  Coating  (NJ,  USA: 

Solec  Inc.),  http://www.solec.org. 


2.  M.  Telkes,  J.  Appl.  Phys.  25,  765  (1954). 

3.  N.  Vatcharasanthien,  J.  Hirunlabh,  J.  Khedari,  and  M. 
Daguenet,  Int.  J.  Sustain.  Energy  24,  115  (2005). 

4.  S.A.  Omer  and  D.G.  Infield,  Sol.  Energy  Mater.  Sol.  C  53,  67 
(1998). 

5.  H.  Xi,  L.  Luo,  and  G.  Fraisse,  Renew.  Sustain.  Energy  Rev. 
11,923  (2005). 

6.  J.  Chen,  J.  Appl.  Phys.  79,  5  (1996). 

7.  G.  Zeng,  J.  Bahk,  J.E.  Bowers,  J.M.O.  Zide,  A.C.  Gossard, 
Z.  Bian,  R.  Singh,  A.  Shakouri,  W.  Kim,  S.L.  Singer,  and 
A.  Majumdar,  Appl.  Phys.  Lett.  91,  263510  (2007). 

8.  A.  Kosuga,  K.  Kurosaki,  H.  Muta,  and  S.  Yamanaka, 
J.  Alloys  Compel.  416,  218  (2006). 

9.  TG  12-4  (TX,  USA:  Marlow  Inc.),  http://www.marlow.com. 

10.  Thermonamic  Electronics  (Xiamen)  Co.  Ltd.,  Xiamen  Tai- 
huaxing  Trading  Co.,  Ltd  (Marketing  Office),  West  Unit  309, 
Guang  Ye  Building,  Torch  Hi-Tech  Industrial  Development 
Zone,  Xiamen,  China  361006. 

11.  Solarbuzz  LLC,  USA.  http://www.solarbuzz.com. 


